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1. Introduction 

In this article we review, in a very elementary way, some of the motivations which are at the 
basis of a class of quantum field theories which have been proposed to describe quantum gravity 
[1, 2, 3, 4]. The base manifold where these models ai^e defined is represented by D copies of the 
Lorentz group, where D is the physical space-time dimension, and the Lorentz group is replaced 
by rotations when the signature of the metric is chosen to be Riemannian. For this reason they 
ai^e known as Group Field Theories (GFT). They can also be viewed as higher rank tensor field 
theories, thus generalising matrix models which describe gravity in two dimensions, and sharing 
some features with noncommutative quantum field theories. The main motivation for us to investi- 
gate such models is that they provide a background independent formulation of quantum gravity in 
which one sums both over topologies and geometries. Indeed, each Feynman graph of a D dimen- 
sional GFT can be dually associated with a specific space-time triangulation. Moreover, it can be 
seen to reproduce exactly spin- foam amplitudes [5], although in a different basis. 

The simplest group field theories correspond to quantization of the BF models, hence to topo- 
logical versions of gravity. The first non-topological model is the Barrett-Crane model [6]. Re- 
cently new spin foam models have been proposed for the quantization of the full 3-i- 1 dimensional 
gravity [7, 8]. These models stem from an improved analysis of the Plebanski simplicity con- 
straints. These new theories could be called dynamical since their propagators, combining two non 
commuting projectors, have non-trivial spectrum. In this paper, after an introductory part dedi- 
cated to motivations and to the discretization of gravity along the lines of lattice gauge theories, 
we focus on such new models, derive the Feynmann rules and outline the computation of quantum 
corrections in the coherent states basis. A detailed account of the results may be found in [9]. 

The paper is organized as follows. In section 2 we review the discrete path integral derivation 
of 2-1-1 gravity in the BF description. In section 3 we introduce the main tools of GFT, which we 
apply to SO{D) BF theories is section 4. We compute Feynmann amplitudes of such models and 
we find that they coincide, for D = 3, with the partition function of 2-1-1 gravity that we have found 
in section 2. This serves as a motivation for the description of quantum gravity in terms of GFT 
models. Therefore we review in section 5 the first order formalism for the gravity action in 3-1-1 
dimensions in the Hoist formulation and we study, in section 6, a GFT model for such an action. 
Finally, in section 7, we describe the stationary phase method used in [9] to evaluate the asymptotic 
large spin regime of certain graphs. 

2. Discrete Gravity in 2+1 dimensions 

In this section we briefly review the first order formalism for gravity in 2 -i- 1 dimensions and 
its discretization, to arrive at a discretized path integral. The partition function so obtained will be 
the basis to introduce and motivate GFT models, which reproduce the same result. 

In 2-1-1 dimensions the first order action of gravity has the form of a BF action with gauge 
group G = S0{2, 1) (S0{3) if the metric is Riemannian) 




(2.1) 
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where e (the B field) is a vector valued one form on the space-time manifold M and F the curvature 

of the connection one form, (O, valued in the Lie algebra ^ = so{2, 1) {so{3) in the Riemannian 
case). In the language of fiber bundles (O is the connection one form of the principal G-bundle, 
while e is a one form on the associated vector bundle, it is therefore vector valued. Once we have 
chosen a trivialization, we can identify the fibers of the vector bundle with the Lie algebra, so that 
e, CO and F are all JSf valued 

e = e'^dx^Zi (2.2) 
CO = (o'^dx^Ti (2.3) 
F = Fl^^dx^Adx^'zi (2.4) 

T/ G -Sf. Therefore, Tr(e AF) is the 3-form constructed by taking the wedge product of the differ- 
ential form parts of e and F and using the Killing form < ., . >= — Tr(..) to pair their Lie algebra 
valued parts. The trace is taken in the adjoint representation. 

In order to discretize this action on a lattice, let us consider a triangulation, T, of the space- 
time manifold in terms of tetraedra, with faces /, edges i and vertices v. We associate to the one 
form e the Lie algebra element 

Ei = [e (2.5) 

To get the discrete analogue of the curvature as a quantity which is still associated to the edges we 
need to introduce a dual lattice, . This is defined on associating to each tetraedron a vertex in 
the middle, v*, to each face an edge which intersects it, to each edge a face whose boundary 
^f^, encircles the edge, to each vertex a closed surface which encloses it. The two-dimensional 
subcomplex contained in such a cellular complex, whose building blocks are vertices, edges and 
faces, draws a graph which we indicate with ^, also known in the literature as dual skeleton. 
Therefore we associate to the connection (O the group element 



hit = Fexp / 



CO (2.6) 



and the discretized curvature is defined as the holonomy of the connection along the boundary of 
/*. Dual faces are associated to edges of the direct triangulation, therefore we have: 

Hf,=He= n (2-7) 

The discretized BF action becomes then 

S{Ei,Hi) = Y^T^rEiHi (2.8) 

ieT 

where the trace is to be taken in the fundamental representation of the group. We can define a 
discretized partition function 

Z[T,n]= f YldEe f Yl dhi^exp{iTriEe JJ hgj (2.9) 
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where the Lie group G is either SO {2, 1) or SO (3), depending on the signature of the space-time 
manifold, while ^ is the appropriate Lie algebra. The integral in the Lie algebra can be performed 
and we get 

We will derive this result from a GFT and this will serve as a motivation to investigate GFT models 
as candidates for quantum gravity. Before going into that, let us recall that (2.10) is also the 
starting point to derive the correspondence of discretized gravity in the first order formalism with 
spin foams. All what is needed is the decomposition of the delta function on the group in terms of 
irreducible representations, labelled by the half-integer j 

5{h)=Y,{2j + l)TrjRJ{h) (2.11) 

j 

where Try is a shorthand for Trvj, the trace in the finite dimensional vector space Vj, base of the 
representation. Replacing into (2.10) and decomposing the representation of the product of group 
elements into a product of representations we finally get 

Z[T*] = nL(2^/* + l)n{67} = A<y (2.12) 

/* jf* 

where all group integrations have been performed. {6j} are the Racah-Wigner 6j symbols. See for 
example [1] for normalization conventions. As indicated, this is exactly the spin-foam amplitude 
Acf of the (dual) graph It represents then the transition from a given space geometry to another, 
as a sum over histories connecting the two. 

3. Group field theory 

In this section we introduce the basic tools of Group Field Theory (GFT) in any dimensions. 
GFT has been first introduced in [1, 2], then further studied by [3, 4], to describe discretized 
models of Quantum Gravity in terms of fields living on the dual lattice of the discretization. We 
first introduce the main ingredients of the theory, such as the notion of field on the group manifold, 
the propagator and the interaction and then we motivate our definitions showing that the quantum 
field theory so defined generates, in a different basis, the spin foam amplitudes of gravity in 2-1-1 
dimensions. 

In our approach all the dynamical content of the models considered is encoded in the propaga- 
tor, while the interaction vertex is fixed by the request that it be the simplest generalization of local 
vertices of quantum field theory. 

We work with Riemannian space-times, therefore the relevant Lie group will be SO{D). In 
three dimensions we will work with its double covering SU (2). In four dimensions we will use the 
splitting S0{4) ~ SU (2) x SU (2). In GFT the field ai^guments live on products of Lie groups 

^ : {gi,...gD) (^[SO{D)f ^ ^{gi,...gD) (3.1) 

with D the spacetime dimension. The propagator is an Hermitean map C : (j) ^C(j) with Hermitian 
kernel C{gi, . . . ,gD;g[, ■ ■ ■ ,g'D)- 

[C(l)]{gi,...,gD) = J dg'i . . .dg'sCigi, . . . ,gD;g'i, . ■ ■ ,g'D)(t>ig[, ■ ■ ■ ,g'D)- (3-2) 



4 



GFT models 



Patrizia Vitale 




Figure 1: A stranded propagator in D=4 with particular orientation; two strands have rjif — \ and the other 
two have rj(f = — 1. 




10 2 6 9 9 8 3 7 



Figure 2: A simple vertex of a 4— dimensional GFT. We have chosen here a particular matching and orien- 
tation for each of the strands. 

It is therefore represented by a stranded line with D strands (fig. 1). The precise form of C is what 
distinguishes the different models. It encodes the dynamical content of the theory. Vertices instead 
only depend on space-time dimensions and are the same for all models. We define a simple vertex 
joining 2p strands in terms of its kernel: this is a product of p delta functions matching strand 
arguments, so that each delta function joins two strands in two different lines. The usual vertex for 
D— dimensional GFT is a 0^+' simple vertex (see fig 2). For instance the SU (2) BF vertex in 3 
dimensions is (with p = 6) 

5int[0] = ^ / ^H'^'^'^ (I>{8i,82,g3)(l>{g4,g5,g6)(l>{g7,gs,g9)(l>{gl0,gn,gl2) K{gi,..,gn), 

with 

K{gu..gn) = 5{g,g^')5{g2g,')5{gegj')5{g,g^^)5{g5gu')S{gig^^) (3.3) 

A stranded graph (see fig. 3) is called regular if it has no tadpoles (hence any line £ joins two 
distinct vertices) and no tadfaces (hence each face / goes at most once through any line of the 
graph). Although final results do not depend on them, it is convenient to introduce orientations in 
terms of incidence matrices: 

• the ordinary incidence matrix e^y which has value 1 (—1) if the edge i enters (exits) the 
vertex v, otherwise. 

• the incidence matrix T] y£ which has value + 1 if the face / goes through the edge i in the same 
direction, —1 if the face / goes through the edge i in the opposite direction, otherwise. 

4. GFT for SO{D) BF theory 

In three dimensions gravity is well described by a BF theory, with gauge group SU (2). The discrete 
partition function has been obtained in (2.10). Such a result is indeed valid for BF in any dimension. 
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Figure 3: The "self-energy" graph '^2 in 4 dimensions, quantum correction to the propagator 

We wish to reproduce such a result with a GFT model. In any dimension the propagator in direct 
space (that is in group variables) is just the projection on gauge invariant fields, 

P(0)=/ dh(^{g,h,...,gDh), (4.1) 

JS0{D) 

= P SO that the only eigenvalues are and 1 (which means that the BF theory has no dynamics). 
P is Hermitian with kernel 

F{g,,...,gD;g\,...g'D)= [ dhfldigMg'ir')- (4-2) 

!=1 

Let us look at amplitudes in direct space. We choose an arbitrary orientation of the lines and faces 
of a graph ^ (which for simplicity has no external legs). Combining the vertex and the propagator, 
integration over all strand variables g,g' leads to 




(4.3) 



where L^, Fg? are the set of lines and faces of ^, respectively. The oriented product fiiefh^''^ 
means that the product of the variables has to be taken in the cyclic ordering corresponding to 
the face orientation (starting anywhere on the cycle). As announced, the amplitude (4.3) neither 
depends on the arbitrary orientation of the lines, nor on those of the faces. This result is valid for 
any D. In order to compare it with our previous results in discrete gravity we have to understand 
the graph 5f as a graph on the dual discretization of space-time. Then we can check that this is 
exactly the result we have found for the partition function of discretized BF (2.10).' Therefore 
the Feynmann amplitudes of GFT with propagator and vertex given respectively by (4.2) and (3.3) 
reproduce exactly the partition function of BF models, which also means, in a different basis, spin 
foam transition amplitudes. 

As we will see in the next section, in 3-1-1 dimensions gravity is described by a constrained 
BF theory. The difficulty in implementing the constraints at the discrete level translates into an 
analogous difficulty in proposing an unambiguous GFT model. 

5. Gravity in 3+1 dimensions: the Palatini-Holst action 

It is very well known that, as for 2-1-1 dimensions, we can describe General Relativity with a 
first order formalism which is, at least classically, equivalent to the second order formalism of the 

'We will understand all GFT graphs as living in the dual triangluation from now on. 
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Einstein-Hilbert action. Differently from the 2+1-D case, it is not just a BF theory, but we have to 
implement constraints. The first order action reads 

S=^l SABCue^ A A [CO] (5.1) 

OTTO Jm 

where ft) = (O^^dx^TAB is the connection one-form, valued in the Lie algebra of S0{3, 1), with 
generators Tab, A = 1,..4 (4x4 antisymmetric matrices). F = F^^dx^ Adx^tAs is the curva- 
ture two-form; it is Lie algebra valued. The tetrads e = e^dx^VA are vector valued one-forms 
such that their wedge product may be regarded as so (3, 1) valued, once a trivialisation chosen 
e Ae = e^eydx^ Adx^TAB- Therefore, similarly to the 2-1-1 dimensional case, the action is entirely 
described in terms of Lie algebra valued forms, the curvature F and the two form e Ae, whose 
wedge product is a 4-form which is integrated on the manifold M; as for their Lie algebra depen- 
dence, on identifying respectively the generators of rotations and boosts with 7a,PQ,fl; = 1, ..,3 we 
can use the nondegenerate bilinear form on the algebra, < . , > i 

< Ja,Pb >l= Sab, < Ja,Jb >l = < Pa,Pb >l=^ (5.2) 

to rewrite the action as 

^1 = ^ / TriC^-AeA/?) (5.3) 
ottG Jm 

with Tri :=< .,. >i. It can be easily checked that it generates the correct equations of motion. 
The peculiar characteristic of the group S0{3, 1) and of its Euclidean version S0{4) is that there is 
another invariant quadratic form on the Lie algebra, Tr2 =< ., . >2, so defined 

<Ja,Jb>2 = <Pa,Pb>2=Sab, <Ja,Pb>2=0 (5.4) 

(actually, eq. (5.4) is the specialization to 3-1-1 dimensions of a quadratic form that would exist for 
any d-i-1). This allows to pair the two-forms F and e A e in a different way so that a new action is 
produced 

S2 = ^ f Tr2(^'A^'AF)= / e"^ Ae'^ AFab[(o] (5.5) 
HkG Jm Jm 

with Tr2 :=< ., . >2- It can be easily shown that this new term can be safely added to the gravity ac- 
tion (5.3) without changing the equations of motion. Although we have no classical consequences, 
we have to expect a different quantum theory. The full action Sh = Si + ^82 is known as the 
Palatini-Holst action [10], while the real pai^ameter 7 is the Barbero-Immirzi parameter [11, 12]. 
Sh can be regai^ded as a constrained BF action 

Sh= f B^^AFAB + ^iB) (5.6) 
Jm 

where the constraint has to implement 

5^fi = £^«e^Ae^ + Va^^. (5.7) 

In the so called Plebanski formulation the constraint is of the form 0(B) °^ <PabcdB^^ AB'^^. The 
action is therefore dependent on three variables, B, (O, <p. Variations w.r.t. to all of them essentially 
give back Einstein equation. 
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The discretization of the 4D BF action S = Jj^^B^^ AFab goes along the same lines as the 3D 
case. We have a triangulation T of the space-time manifold with a 4-D simplicial complex, and a 
dual lattice T* with its 2-dimensional subcomplex which identifies our graphs. The B field now 
being a 2-form, has to be integrated over faces / to produce so{4) (Lie algebra) elements, Bf. The 
discretized curvature is defined as in 3D to be the holonomy of the connection along the boundary 
of /* (2.7). Dual faces are now associated to faces of the direct triangulation, therefore we have: 

SiBf,Hf) = £ TrBfHf (5.8) 

When writing the partition function for such a discretized action, upon integrating over the Lie al- 
gebra, we end up with the same result as in (2.10) with the gauge group now SU (2) x SU (2). When 
using the Peter Weyl theorem to decompose the group delta function on irreducible representations, 
we obtain a spin foam amplitude which is the analogue of (2.12), with 6j symbols now replaced by 
15j symbols (see [2] for normalization conventions) and j = {]+,]-), according to the splitting of 
SO{A). 

To treat the constrained action (5.6) we have to include the Plebanski constraint. There are 
various proposals on how to implement it at the discrete level. We have considered just one of 
them, the EPRL spin-foam model [8], in the Freidel-&asnov formulation [7], which uses coherent 
states. This is the argument of the next section. 

6. GFT approach to the EPRL-FK model 

In our approach all the dynamical content of different models is encoded in the propaga- 
tor, while the vertex remains the same as in BF theory (3.3). Therefore the constraint is di- 
rectly imposed on the propagator (4.2). This is better understood in the coherent states basis. 
In 4 dimensions we use the splitting of S0{4) as the product SU (2) x SU (2) and we decom- 
pose the group elements as g = with g± G SU (2). In 3D we use SU (2) coherent states 
\j,g >= g\jj >= Y,m li)"^ > ^i'-'^ ^^^Hg) the 2j-i-l dimensional in^educible representa- 
tion of g. The identity is decomposed over coherent states as 

^j = ^j dg\j,g>< j,g\=dj dn\j,n><j,n\ (6.1) 

Jsu(2) Jgih=s- 

with \i^n >= R^^\gn)\j,j >■ In 4 dimensions we indicate with j = the eigenvalues of 

the Lie algebra generators J in each SU (2) component and we work with SU (2) x SU (2) coherent 
states > >. The identity is therefore 

l; = d;,d; / I 7'+, ?I+ > Cg"! /- , >< /+, ^ < J- , 1 (6.2) 

As an intermediate step towards the full 4-dimensional case we rewrite the D-dimensional BF 
propagator (4.2) in a different way (D=3, respectively D=4). Since = P we introduce two gauge- 
averaging variables, m,v G G at the ends of the propagator, {u on the side where e,,^ = — 1 and v 
on the side where = +1)- G is either SU (2) or SU (2) x SU (2) depending on the dimension. 
Between these two variables we insert the partition of unity (6.1) (resp. (6.2)) 
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where, to simplify the notation, we omit to indicate the representation expUcitly and use g = R^j^ (g) 
for all group elements, from now on. This is equivalent to (4.2) upon redefining h = v^^u. The 
index / labels the strands of the propagator according to the face to which they belong. In other 
words, on introducing 

I = ®jf Ijf (6.4) 

we have used 

(6.5) 



This redundancy will be useful in a while to better understand the constrained propagator. 

The amplitude of a given graph ^ with a certain number of vertices and propagators may be 
written as a factorization over faces 

f Yl ^"^'^^^ n -^Z (6.6) 

We assume for simplicity that all faces are closed. Then, for each given face / of length p, we 
number the vertices and lines in the (anti)-cyclic order along the face as £i , vi • • • £p,Vp, with ip+i = 
l\. Combining vertices (delta functions for the strand group variables, g,,g-) and propagators, and 
integrating over the strand group variables, gi,g\ we obtain 

where is if £4^,, = +1 and ui^ if e^^,,^ = -1. (Remember that in D=4 dy = dj^^dj^J). 

6.1 The EPRL/FK GFT 

The EPRL/FK propagator has a structure similar to (6.3) but with replacement of I by a non- 
trivial projector which does not commute with P. This implies that it is not possible to recombine 
u and V in a single gauge averaging variable h. It implements in two steps the Plebanski constraints 
with a non-trivial value of the Immirzi parameter 7. We only consider < 7 < 1 where the EPRL 
and FK models coincide. Starting from the expression (6.3) of the BF propagator in the coherent 
states representation, the first step adds the constraint ]+/ j- = (1 + 7)/(l — 7), = The 
second step replaces in each strand of (6.3) the identity \j by a projector tJ whose definition is 



= dj++j- [^y>-/;>+=(i-r)/(i+r)] / MJ+^ri > ®\j-,n >< j+,n\® < j^,n\ 



(6.8) 



Grouping the four strands of a line defines an operator, T^, that acts separately and independently 
on each strand of the propagator: 

T^ = (Bjf^UTl (6.9) 
so that the EPRL/FK propagator, to be compared with (6.5), is 

C = PT^P (6.10) 
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that is 



C{g,g) = / ^^"^/vnE[V/i/+=(i-r)/(i+r) 
■' /=! if 




(6.11) 



Tr;/+®i/- (§/) ^ V ^ 1 7/+ , «/ > 1 7/_ , Uf) {jf+ , «/ 1 ® < 7/_ , 1 j , 

with aj = dj^dj_ , pj = dj_^^j_ . We state, without proof, a number of properties of the propagator, 
for more details see [9, 13] 

• The operator C is hermitian 

• is a projector 

• Since the propagator is hermitian, Feynman amplitudes are independent of the orientations 
of faces and propagators. 

• Since and P do not commute, the propagator C can have non-trivial spectrum (with eigen- 
values between and 1). 

• Slicing the eigenvalues should allow a renormalization group analysis. This is why we call 
this kind of theories dynamic GFT's. 

• Since is a projector, the propagator C of the EPRL/FK theory is bounded in norm by the 
propagator of the BF theory, as well as Feynman amplitudes. 



6.1.1 EPRL/FK Amplitudes 

Combining the propagator and the vertex expressions, the integrations over all g,g' group variables 
can be performed explicitly, leading to the amplitude of any graph This amplitude is given 
by an integral of a product over all faces of the graph as in (6.6), but the amplitudes for faces 
are different. To compute these face amplitudes we distinguish between closed faces (no external 
edges) and open faces (which end on external edges). For example, the self-energy graph in fig. 3 
has 6 closed faces and 4 open faces. Using the same numbering of the p edges and vertices along 
a closed face, its amplitude is given by 

^f= j f\{dScM:)I.<^j,J^i^.+®i^Ai''t^^^^^^^ (6-12) 

where the constraint on j+,j- is implicitly understood from now on and Vy are, the vertices. To 
perform the integrals over the strand variables g,g' we use the standard results 

i5(7,/)5-5,f (6.13) 

^d{j,j')s'"Pe„, (6.14) 
"i 



J dgR^j\,'"{g)R'^j'^P,{g) = 
J dgR^J>\{g)R'^j>,{g) = 
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with R^j^'"„{g) the matrix element of g in the 2j + I dimensional representation and R n'"{g) = 
R^^''"'n{g^^) which imply 

j dgTvjAgTvj,Bg-' = j5{iJ')TvjAB (6.15) 

j dgTTjAgTvj,Bg = j5UJ')TrjAeB^e^ (6.16) 

and £ G SU{2), 




(6.17) 



Thus we arrive at 

= L^jf^^jf.-jf-UL (^jH:::iTj) , (6.18) 
jf 

with £p^i = l\. It is interesting to notice that we recover the SU (2) BF theory in the limit 7 — )• 1 
(see [9] for details). For open faces (which end on external edges) we obtain a slightly modified 
expression but we also find that their contribution is in^elevant in the large spin approximation. It 
can be shown that the graph amplitude obtained by replacing (6. 18) in (6.6) is the same as the spin 
foam amplitude found in [8] when passing to spin representation. 



7. The stationary phase method 

In [9] we have evaluated the degree of divergence of the graphs for GFT models using the 
stationary phase method. Here we give a brief account of the calculation, which is however quite 
technical, and we state the results. We refer to the original paper for more details. To this analysis, 
we have to rewrite the graph amplitude (6.6), with face amplitudes depending on the model, in the 
form 

f U'^hlldnexpiY^jfSfiKn]}, (7.1) 

;><A •' f 

with ^ = ■J^{j) and A some ultraspin cutoff. To evaluate the superficial power counting, we set 
if = jkf with kf G [0, 1] and use the stationary phase method to derive the large j behavior of the 
summands. This calculation has been performed both for BF models and for the EPRL/FK model. 
If the action is complex but has a negative real part, the contribution to this integral are quadratic 
fluctuations around zeroes of the real part of S which are stationary points of its imaginary part, 
otherwise the integral is exponentially suppressed as j — )• 00. 

7.1 The saddle point method for 3D BF models 

For 3D BF models the degree of divergence of the graphs has already been calculated in the 
literature [15, 16]. We re-derive such results to test our method. To write the action Sf we introduce 
the projector 

\n){n\ = ^{\ + a-n), (7.2) 

with \n) = \n,^), and observe that 

{n,j\g\n',j) = {n\g\n'fJ (7.3) 
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so that 



Sf[h,n]=kf\ogTr {W hf^) {\ + o ■ n) 

e.edf 



(7.4) 



Since the action is the logarithm of the trace of the product of a unitary element and a projector, its 
real part is negative and maximal at /z^ = 1 . To perform the saddle point expansion, we expand the 
group element to second order as 



^2 

h = l-^+iAro + 0{Aj) (7.5) 



with A € su{2). Also 



nf = nf+^f-^nf^ + 0{^}), with •<§/ = 0. (7.6) 

(because n^^ = 1 up to third order terms). Let us consider a face with edges £i,...,£p, then to second 
order 

Hh^/'^ = l-^+iG-Af-iG-^f (7.7) 

with 

L and <!>/= £ 7]e,,file,,fAe, ^Ae,. (7.8) 

l<a<p l<a<b<p 

this implies 

SflAiAf] = 2kf{inf ■Af-^ + i^?^^ +i^f.Af + inf ^f] (7.9) 
and we have to estimate 

J YldAelld^f exp2;X^/{i«/ • A/ - ^ + ^^^l^ .Af + iuf ■ ^f] (7.10) 
as j — > oo. 

We do not integrate over the vectors n/. They have to be chosen so that they are extrema of 
the imaginary part of 5. 

Because all terms except the first one Y^f^f^f 'Af are of second order, the imaginary part is 
stationary if and only if 

Y,ikfnfAf = Y,ini,fkfnfAi = VA^ e M^ (7.11) 

which amounts to the closure condition 

Y,ne,fkfnf = 0, ye (7.12) 

/ 

This is the requirement that, in the semi-classical limit, the vectors jfUf are the sides of a triangle 
(resp. the area bivectors of a tetrahedron) that propagates along linD = 3 (resp. D = 4). 
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The solutions of the closure condition (7.12) range from non degenerate to maximally de- 
generate. In three dimensional (resp. four dimensional ) BF theory, a solution is said to be non 
degenerate if all the tetrahedra (resp. 4-simplices) corresponding to the vertices of the graph have 
maximal dimension. At the opposite end, we say that a solution is maximally degenerate if all the 
vectors «y are proportional to a single one n^, 

nf = Ofno with a/G{-l,+l}. (7.13) 

In both cases we find 

i/^~A^^"^'' (7.14) 

with r the rank of the Lx F incidence matrix rjtf. 

In particular for the self-energy graph in fig. (3) we find 

i%~A'^ (7.15) 

with F = 6, the number of closed faces. We do not include open faces because they do not con- 
tribute to the divergence of the graph, as we will argue below. 

7.2 The self-energy in the EPRL/FK model 

Let us come to the EPRL/FK model. An analysis for generic graphs is possible although 
technically difficult. Here we only consider the self-energy graph ^2 of fig. (3), which represents 
the first quantum correction to the propagator. Moreover, we limit the analysis to non-degenerate 
configurations, as in [14]. For completely degenerate configurations the stationary phase method 
only provides upper bounds. The issue is discussed in [9]. 

The self-energy graph '^2 has 4 open faces and 6 closed faces with two edges. We label the 
internal propagators with an index a ranging from 1 to 4 and orient them in the same direction. We 
label the 6 closed faces with pairs of indices (a, ft), a <b. We do not worry about open faces in our 
analysis, which technically amounts to put to zero the external spins. However, the result should 
also hold for finite non-zero spins on the external faces. Indeed, since the latter remain finite, the 
contribution of the external faces to the action can be neglected as j — )• 00. The amplitude of the 
self-energy graph reads then 

a a<h 

where, from (6.18) the face amplitude reads 

Kb = Y.^udj_pj / dnabdn'^,^{j+nah\ua+u^l\j+n'^h){j+K<b\'^b+Vai\j+nab) 
j 

X {i-nab\ua-u^l_\j-n^^,,){i-n[,i,\vh-v^l\i-nab). {1 Al) 

In order to perform a stationary phase analysis we rewrite the graph amplitude as 

^^2=L[ n^"^ n^^^^ lldnill{{dj^fdj.dj-exp{jS}+jSj}} (7.18) 

jfJ a a r f ^ ^ 
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with 7j = jy^kf, kf G [0, 1] and j large. The face action for f = abis then 

Sf = 2r%log{{nf,,\u^{u^)-'\nfjMnfMvU^ar'M} (7-19) 
We employ the saddle point technique around the identity 

uf = l-^^+ia-Af + 0{A^)\v^ = l-^^+ia-Bf + 0{B^f (7.20) 
and use the result (7.2) so that the action at the identity for the face f = ab reads 

S±[l,l,«,] = 7±^„,log{^^i"^^} (7.21) 

which is negative except for n/^ ^, = Ufj, = nf. 

We perform the expansion of the coherent state around a unit vector common to all the strands 
of the face 

ni = nf + ^i-^^nf + 0{^i)\ with «/-^,- = 0, (7.22) 

otherwise the integral is exponentially damped. It is convenient to perform the following change 
of variables 

A^=Aa±r^Xa and B^=Ba±r^Ya (7.23) 

Terms linear in A* and only involve A and B, while in the quadratic terms, the pair of variables 
A and B on one side and the pair X and Y on the other side decouple. 

Thus we observe that, at the level of the quadratic approximation, we can separate the action 
into a SU(2) BF action (variables A and B) and an ultralocal potential that only involves uncoupled 
variables attached to the vertices (variables X and Y). Moreover it can be shown that this result is 
true for a generic graph. 

We perform the Gaussian integration over the two dimensional vector Xf = ^f,a — ^f,b, and 
obtain 



=2% = 

Jf 



L { / n^^« n^^« lld^fexpjSBF{A,B,^) (7.24) 
j Y\dXatxpiQ{X) X / Y\dYatxpjQ{Y)\ (7.25) 



with i,f = (^y a + ^j\h- The BF-like action is 

5bf[A,S,^] = £;:,i{-^[«/A(A„-Afe + Sfc-S„)]^ (7.26) 

a<b 

+inab ■ {Aa -Ah+Bh-Ba) +inab- {Aa AAb+BbABa) (7.27) 
+i U • {Aa -Ai, + Bh - B,) } (7.28) 

while the ultra local terms are 

Q[X] = Y+r L kab{ [nab A iXa-Xh)f + inab ■ (X„ AXi,)}. (7.29) 

a<b ^ ^ 
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The Gaussian integral over tlie variables A and B can be evaluated using the same techniques as for 
the pure BF models, thus yielding 

IX{dAaY{dBaY\d^ftx^jSBF{A,B,^) ~ (7.30) 

■J a a f 

as j — )• oo. As for the Gaussian integrals over the independent variables Xa and Ya, we have 

lldXaexpjQiX) ~ j-— (7.31) 

a 

with Q[X] the quadratic form in Eq. (7.29) and a similar expression for Q[Y]. The rank may be 
computed to be rank(Q)= 9 so that the Gaussian integral over X,Y yields a power of j'^^/^ each. 
Therefore, we obtain the power counting for the self-energy with non degenerate configurations as 
follows 

E /'xr'xr'x(r^/2)^~A^ (7.32) 

6 independent spins ~ j <A 

where 

• y^"* arises from a dj+dj- ~ for each of the 6 faces 

• dj ~ i for each of the two strands in each face 

• j^^ results from the Gaussian integration over the 6 variables Xf = i^f.a — ^/m) 

• from the integration over A and B 

This reproduces the result of [PRS], with non degenerate configurations. 
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